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/\ ' Abstract 

We analyze in detail the highly damped quasinormal modes of d-dimensional 
Reissner-Nordstrom black holes with small charge, paying particular attention to the 
large but finite damping limit in which the Schwarzschild results should be valid. 
In the infinite damping limit, we confirm using different methods the results ob- 
tained previously in the literature for higher dimensional Reissner-Nordstrom black 
holes. Using a combination of analytic and numerical techniques we also calculate 
the transition of the real part of the quasinormal mode frequency from the Reissner- 
Nordstrom value for very large damping to the Schwarzschild value of ln(3)rf,/i for 
intermediate damping. The real frequency does not interpolate smoothly between 
the two values. Instead there is a critical value of the damping at which the topology 
of the Stokes/anti-Stokes lines change, and the real part of the quasinormal mode 
frequency dips to zero. 



1. Introduction 



There exist by now many calculations of the highly damped quasinormal modes (QNMs) 
of a large variety of black holes. Since the highly damped modes are not directly observ- 
able, these calculations are mainly motivated by the proposal, originating with Hod Ij and 
more recently propounded by Dreyer[2] in the context of loop quantum gravity, that these 
modes may be providing information about the semi-classical quantum spectrum of black 
holes. This proposal was prompted in part by the apparent universality of the real part 
of the frequency in the highly damped limit, as well as by its special numerical value for 
Schwarzschild black holes, namely: 

un^--^ \n{3)nj,/h, (1.1) 

where Tbh is the Hawking temperature of the black hole. The frequency and the damping are 
proportional to the temperature on dimensional grounds for single horizon black holes, but 
the coefficient ln(3) allows, on application of the Bohr correspondence principle, an elegant 
statistical interpretation U 0] of the resulting Bekenstein-Hawking entropy spectrum. We 
refer the reader to the literature for details|3]. Although this special value has been shown 
to be valid for a large class of single horizon black holesj^linilZllH] this simple relationship 
between frequency and temperature breaks down for multi-horizon black holes. 

It was shown first by Motl and Neitzke[5j and then confirmed by Andersson and Howls jH] 
that for Reissner-Nordstrom (R-N) black holes the real part of the highly damped QNM 
frequency in four dimension approaches ln{5)Tbh as the charge goes to zero. The general 
validity of this result was verified by Natario and SchiappajTU], who calculated the R-N 
highly damped spectrum in arbitrary spacetime dimensions. The apparent contradiction 
between the zero charge limit of the R-N case, and the Schwarzschild value was explained 
heuristically by the authors of [S]. They noted, that, while the ln(5) represented the correct 
R-N result for very high damping, one expects an intermediate range of damping for which 
the Schwarzschild value of the real frequency is correct. Order of magnitude arguments 
suggest this range in 4-d to be: 

1 « IcupM^ « MVgS (1.2) 

The purpose of the present paper is two- fold. First, we calculate the highly damped 
QNMs of d-dimensional R-N black holes with the small charge using the techniques of 
Andersson and Howls 9], thereby confirming the results of Natario and Schiappa JU], who 
used the monodromy method of Ref . jHJ [TT] . Secondly, and more significantly, we use a 
combination of analytic and numerical techniques to analyze the limit of large but finite 
damping where the Schwarzschild limit is approached. For four and five spacetime dimen- 
sion R-N black holes we explicitly calculate the spectrum in the transition region from 
ln(3 + 4cos(^^7r))T6/i for very large damping to the Schwarzschild value of \n(3)Thh- The 
real frequency does not interpolate smoothly between the two values. Instead there is 
generically a critical value of the damping at which the Stoke/anti-Stokes lines change 
topology and the real part of the frequency dips to zero. This behaviour seems to be the 
analogue of the algebraically special frequencies (for a recent discussion see Berti^SJ that 
mark the onset of the high damping regime. We provide general arguments that suggest 
the persistence of this behaviour for all spacetime dimensions. 



The paper is organized as follows. In Section 2 we describe the general formalism. 
Section 3 summarizes the results of the calculation in the large damping limit for metric 
perturbations (tensor, vector, and scalar) of Schwarzschild and R-N. The former reviews 
well known results while the latter provides independent confirmation of JU]- Section 4 
presents the details of the calculation for axial perturbations in the intermediate region 
in four dimensions, while Section 5 does the same for five dimensions. (Tensor and scalar 
perturbations cannot be done using our numerical techniques, because a crucial term in 
the effective potential vanishes precisely in that limit.) The final section describes general 
arguments for the persistence of the qualitative features of our results in higher dimensions, 
and closes with some conclusions. 



2. General Formalism 

It was shown by Ishibashi and Kodama in ^^1, ^M^ ^^^ ff3 that various classes of 
non-rotating black hole metric perturbations in a spacetime with dimension d > 3 are 
governed generically by a Schrodinger wave-like equation of the form: 



d'i: 



+ 



uj^-V{r) ^ = , (2.1) 



for perturbations that depend on time as e"*"^*. The Tortoise coordinate z is defined by 

dr 
d. ^ ^ . (2.2) 

where /(r) is related to the spacetime geometry, and is given by 



/W = l-;^ + -^-Ar^ (2.3) 



The ADM mass, M, of the black hole is related to the parameter /i by 



SnGd 
where An is the area of a unit n-sphere. 



M = <''-y^-V , (2.4) 
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The electric charge, q, of the black hole is, 

(rf_2)(d-3)^2 



^" = / ^ix • (2-5) 

■p / n+l ^ 



q' = '—-^ ^-e' , (2.6) 



while the value of the cosmological constant. A, is given by 



(d-l)(d-2)^ , , 

A = ^ '-^ ^A . (2.7) 



In this paper we only deal with asymptotically flat black holes where A = A = 0. The 
units we will be using are c = Gd = I- 

The effective potential V{r), was found explicitly by Ishibashi and Kodama fini IT^ 113] 
for scalar (reducing to polar at c/ = 4), vector (reducing to axial aX d = 4), and tensor 
(non-existing at d = 4) perturbations. The effective potential for tensor perturbations has 
been shown [TBI HTj to be equivalent to that of the decay of test scalar field in a black hole 
background. The effective potential is zero at both the horizon {z -^ — oo) and spatial 
infinity {z -^ oo). In the case of QNMs, the asymptotic behavior of the solutions is chosen 
to be 

nz) ~ I ^+^uz as 2 ^ 00 (x ^ 00) , ^^■^' 

which represents an outgoing wave at infinity and an ingoing wave at the horizon. Since the 
tortoise coordinate is multi-valued, it is more convenient to work in the complex r-plane. 
After rescahng the wavefunction ip = ^/v7 we obtain 

+ R{r)^ = 0, (2.9) 



dr'^ 



where 



with 



u' 



Rir) = y^ - Uir) , (2.10) 



V{r) If" Iff 
P 2 f 4 

Here prime denotes differentiation with respect to r. 
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-«-7^-^T-i 7 ■ (^■"' 



3. Large Damping Limit and Universality 

We now consider the QNMs in the infinite damping limit where 

\uj^\ ^ [Im cuf ^ 00 . (3.1) 

Since we will be using complex analytic techniques, in principle the behaviour of U{r) on 
the entire complex plane may be relevant. However, in the infinite limit case, the io^ j P 
term in K{t^ will dominate f/(r) everywhere, unless one of the terms in f/(r) diverges. This 
can only happen at the origin, or at the horizon. However, since oj"^ / P also diverges at the 
horizon, it will dominate there as well in the large damping limit. Thus, only the dominant 
term of U{r) near the origin is relevant in this limit and R{r) can be approximated on the 
entire complex plane by 

where J is determined by the type of metric perturbation. In particular, for different 



perturbations we find: 



J' -I 



— 1 Tensor Perturbation (0 = and n > 2) , 

n{n — 2) Tensor Perturbation (0 7^ and n > 2) 

4n^ — 1 Vector Perturbation (6' = 0) , 

3n(3?2 — 2) Vector Perturbation (6* 7^ 0) , 

—1 Scalar Perturbation (0 = 0), 

ra(ra — 2) Scalar Perturbations (0 7^ 0) , 



(3.3) 



where n = d — 2. The universality of the infinite damping QNM frequency can, to a large 
extent, be attributed to the simple form of ()3.2|1 in this limit. 

In the WKB analysis, the two solutions to Eq. ()2.9j) are approximated by 



^W(x) = Q-^/\x) exp [+t /," Q{x')dx'] 
Q-^/^{x) exp[~i 1^"= Q{x')dx'] 



*?(a; 



where 



Q'ir) = R{r) - ^ - ^ 



4^2 



(3.4) 



(3.5) 



is shifted by l/(4r^) in order to guarantee the correct behaviour of the WKB solution at 
the origin. In Eq. ()3.4|1 . t is the simple zero of Q^. 

Given the function Q, it is possible to determine the WKB condition for asymptotic 
QNMs. The methodology that we adopt has been explained in details in [Hj. First, one 
determines the zeros and poles of the function Q and consequently the behavior of the 
Stokes and anti-Stokes lines in the complex r-plane. Stokes lines are the lines on which the 
WKB phase (/ Qdr) is purely imaginary and anti-Stokes lines are the lines on which the 
WKB phase is purely real. It turns out that for Schwarzschild black holes in any spacetime 
dimension, there are two unbounded anti-Stokes lines which extend to infinity on either 
side of a bounded anti-Stokes line that encircles the event horizon. In the case of R-N 
black hole, which has two horizons, the unbounded anti-Stokes lines are again on either 
side of a bounded anti-Stokes line encircling both the event horizon and the inner Cauchy 
horizon. Inside this loop a second bounded anti-Stokes line encircles the inner Cauchy 
horizon. The structure of anti-Stokes lines for Schwarzschild and R-N black holes in four 
and higher dimensions can be found in pUj . 

To determine the WKB condition on the QNM frequency we start on an unbounded 
anti-Stokes line where the solution is known due to the boundary condition at infinity: 



'^initial — ^c 



(3.6) 



It is relatively simple to determine how the solutions change as we move along an anti- 
Stokes line. Moreover, there are rules that specify how the solution changes when one 
moves from one anti-Stokes line to another while crossing Stokes lines. The procedure is to 
move along a contour that goes from the known solution at infinity along the unbounded 
anti-Stokes line, changes to the bounded anti-Stokes line(s), loops around the event horizon 
and finally returns to infinity along the same unbounded anti-Stokes line that it started 



from. This yields a final wave solution ip final that is related to the one we started from as 
follows: 

i^ final = X"^ initial , (3.7) 

where x is the monodromy of this contour. Since this contour circles only one pole, namely 
the horizon, the monodromy of this contour must be the same as the monodromy of 
a contour in the vicinity of the pole at the event horizon. This latter monodromy is 
determined by the boundary condition at the horizon, so that we are left with a consistency 
condition that determines asymptotic QNMs. 

Applying the above method to Schwarzschild black holes in four and higher dimensions 
one arrives at the condition 



^2ir 



— 1 — 2 cos ( —IT 



(S.t 



where the appropriate value of J is given in ()3.3|) . F is the integral of the function Q along 
a contour encircling the pole at the event horizon in the negative direction. Using the 
residue theorem to evaluate T yields 



r = ® Qdy = —2'Ki Res Q = —2Txi 

J r=\rh\ 



1 



n — 1 



(2^) 



l/(n-l) 



UJ 



(3.9) 



We can now combine Eqs. (j3.3|l . (|3.8j) . and ()3.9|1 to evaluate a;: 

47r 1 

■ {2 jj,)"^^ UJ = \n{3) + {2k +l)TTi as k -^ oo 



n 



(3.10) 



This result is valid for all types of metric perturbations in four and higher dimensions. 

In the R-N case where 6' 7^ 0, the same procedure yields the following general WKB 
condition for four and higher dimensions: 



,2ir. 



1 



1 + cos 



J 



2n 



-TT 



:i + e- 



-2iri 



(3.11) 



where Fe and Fj are the integral of the function Q along a contour encircling the pole at 
the event horizon and the inner Cauchy horizon respectively in the negative direction. We 
can also evaluate these integrals using the residue theorem, and the results are 



-27Ti Res Q = —ni- 



and 



Fj = —27Ti Res Q = iri- 

r=|r_| 



[n-l)y/JF^^ 



r^u 



(3.12) 



(3.13) 



Here r+ and r_ are the locations of the event horizon and the inner Cauchy horizon 
respectively, which are given by 



r±=U±V/i2-^2 



(3.14) 



Since in this paper we are interested in the Schwarzschild hmit of Eq. ()3.17|) . let us find uj 
in the hmit where ^ — !► 0, we note that in this hmit: 

- 27r? , , 1 , , 

-c<;(2/i)— , (3.15) 



n 
and 



n — 1 
Equation ()3.11|) then takes the simple form 




(3.16) 



exp ( — —i2^i)^^uj\ ^ -3 - 4 cos ( vr) . (3.17) 

\n — 1 ) \ln — I J 

The constant J for tensor and scalar perturbations is ri — 1, while for vector perturbations 
it is 2n. Even though J is different in vector perturbations, but for these perturbations we 
can write 

2n-l 2n-l 2n - 1 ' ^ ^ 

which produces the same result in Eq. ()3.17|) as tensor and scalar perturbations for which 

^ "-' (3.19) 



2n-l 2n-l 
Thus, in the small charge limit, the R-N black hole QNM frequency is: 

47r 



(2u)"-ico' ^ m 3 + 4cos n ] 

n-V ^^ [ \2n-l J 



+ {2k + 1)711 as fc -^ oo . (3.20) 



Clearly this does not correspond with the Schwarzschild result in Eq. ()3.10|) . The difference, 
as pointed out in |9j is due to the order in which the limits |a;| ^ cxd and ^ — > are taken. 
This issue will be addressed in detail in the following sections. 

Note that the results in this section are consistent with the results obtained by Natario 
and Schiappa in 10|, where the reader can find a detailed calculation of the highly damped 
QNMs based on the different methodology first used in j3] and [TT| . 



4. Schwarzschild Limit of Reissner-Nordstrom Black Holes 

It was shown in the previous section that the infinitely damped QNM frequencies of 
R-N black holes do not approach the Schwarzschild result in the limit when the black 
hole charge goes to zero. This is a potentially puzzling issue that deserves more careful 
investigation. The resolution was suggested in |9|: there should be an intermediate regime 
of large but finite damping in which QNMs resemble the Schwarzschild result. We now 
investigate this issue using a combination of analytic and numerical methods to carefully 
track the modification in the structure of the Stokes and anti-Stokes lines in the complex 
plane as we approach the Schwarzschild limit. We will of necessity concentrate only on 
vector perturbations because in tensor and scalar perturbations the constant J is zero for 
Schwarzschild black holes where ^ = 0, which makes the numerical analysis problematic 
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Figure 1: The zeros of the function Q for pure gravitational axial perturbations in four 
spacetime dimensions for different values of |cc;|. 



in the Schwarzschild limit. Note that the analytic calculations stay away from J = until 
the end, and then take the limit J — i> in the final expression for the frequencies. 

In the case of R-N black hole, for vector electromagnetic-gravitational perturbations 
one has two distinct effective potentials 



K± = /(r) 
where 

Setting A = we get 
1 f 



n{n-2) 4/(/ + n-l)+n(n-2) /?±/i n{5n - 2)9'^ 



4 

n^ + 2 



47-2 



r.n+1 



Aj,2n 



± W(n2 - 1)2 + 2n{n - !)[/(/ + n - 1) - n] 



9^ 



(4.1) 



(4.2) 



„ ,, . , 2 Al(l + n-l)+n(n-2) Al(l + n - 1) + 3n^ - An ~ (3± 
Rv±ir) = p[^ ^2 + ^^1 /^ 

n^-l-2(3± 2 n(7n-5) + f(n-2) + /(/ + n-l) 



r2n 



-fi - 



^2n 



3n{n ~ 1) + mn - 2) - (3^ ^ 



2,.,. ^j ^±_^^2 ri(n-l) + f(5n-2)^,] ^ ^^^^^ 



r.3n— 1 



„4n-2 



Since we are interested in the Schwarzschild limit, throughout this paper we take 9/fi ^ 1, 
where the coupled electromagnetic and gravitational vector perturbations approach pure 
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Figure 2: Anti-Stokes lines emanating from the two zeros before (a) and after (b) they 
coincide on the real axis. 

electromagnetic and gravitational perturbations. In this limit, for electromagnetic waves 
we have /5+ ~ ra^/2 — 2 and for gravitational waves we have /3_ ~ — 3n^/2. In order to 
do numerical analysis, let us concentrate first on four dimensional spacetime where the 
function R for pure electromagnetic and gravitational perturbations is 

(4.4) 
Here q and M are the charge and mass of the black hole respectively and j is the spin 
of the perturbation which is 1 for electromagnetic and 2 for gravitational perturbations. 
Note that in the infinite damping limit, where — Im uj —>■ oo, the zeros of the function R 
approach the origin of the complex plane and as in Eq. ()3.2|) one can always take 

R{r) - ^ - ^ • (4.5) 

On the other hand, for large but finite damping, where — Im cu ^ Re c^;, the zeros of the 
function R are not infinitely close to the origin of the complex plane. In this case we need 
in principle to consider all the terms in the function R. The terms in Eq. fl4.4|) of order r~^ 
and r~^ inside the curly bracket are in fact negligible as long as the zeros of the function R 
are located at a radius much smaller than 1, but for the sake of accuracy, we nonetheless 
keep all the terms in Eq. ()4.4p in our numerical calculations. 

Fig. ^ shows how the zeros of the function Q^ = R — ^, for pure gravitational axial 
perturbations (j = 2), move in the complex plane as we change the value of \u!\. (In our 
calculations we always assume that |u;| = |Im cjI.) In this figure we have taken 1 = 2, 
2M = 1, and q = 0.01. As |a;| decreases, four of the zeros move away from the origin 
and approach the Schwarzschild limit in which the zeros are located at tq ~ Iroje'-^'^"'"^-''^/^, 
where A; = 0, 1, 2, 3. The remaining two zeros, located on the right side of the origin, 
initially move away from the origin and then they move toward the positive real axis, 
where they eventually coincide at a point to the right of the inner horizon. In Fig. ^ 
the two zeros coincide when the damping is approximately equal to 5.3 x 10^. After they 
coalesce, one zero moves toward the inner horizon and the other one moves away from the 
inner horizon. This motion does not continue indefinitely. In Fig. ^ the two zeros on 
the real axis asymptote to r = 0.00009122 ■ ■ -, which is just inside the inner horizon, and 

9 




Figure 3: Schematic presentation of Stokes (dashed) and Anti-Stokes (sohd) hnes in the 
complex plane before the two zeros coincide. The open circles are the zeros and the filled 
circles are the poles of the function Q^. 

r = 0.00017133 • • •, which is just outside the inner horizon. Thus, two of the zeros stay 
relatively close to the inner horizon while the other zeros move far away from the origin. 
These four zeros finally enter a region in the complex plane far enough away from the origin 
so that the terms with the charge q in the functions /, R, and Q become negligible. This 
happens in the range 



(j2-l) + 3/4 , ,2 4(j 



(2M) 




(4.6) 



For q/M <^ 1, this nonetheless allows many orders of magnitude for |Im u;|. In this range, 
we can approximate Q^ to be 



Q\r) ~ - jo. 



lU' 



r 



3/4] (2M)^ 



1 



(4.7) 



which is exactly identical to what we have for Schwarzschild black holes. 

It is important to look carefully at what happens when two of the zeros merge on the 
real axis. Fig. Efa), shows the two zeros just before they coincide on the real axis. The 
two anti-Stokes lines to the left encircle the inner horizon and the two anti-Stokes lines to 
the right connect to the neighboring zeros on either side of the positive real axis. After the 
zeros coincide the angle between the anti-Stokes lines on the right hand side, which was 60° 
before they coincide, has changed to 120° as one can see in Fig. |21 During this transition, 
the lines must break apart. The anti-Stokes lines to the right detach from their neighboring 
zeros and encircle the pole at the origin. The schematic behavior of the Stokes and anti- 
Stokes lines before and after the two zeros coincide on the real axis are shown in Figs. |21 
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Figure 4: Schematic presentation of Stokes (dashed) and Anti-Stokes (sohd) hnes in the 
complex plane after the two zeros coincide. The open circles are the zeros and the filled 
circles are the poles of the function Q^. 

and E] respectively. These schematic figures are based on numerically generated diagrams 
of Stokes and anti-Stokes lines in the complex plane. An example of such numerically 
generated diagrams is shown in Fig. El which plots the Stokes and anti-Stokes lines for the 
gravitational perturbations (j = 2) with / = 2, 2M = 1, q = 0.01, and |ci;| = 5 x 10^. Note 
that Fig El shows that at this point one of the zeros is located inside the inner horizon. 

The numerical recipe we have written is operating in the following way: The input 
parameters are j, M, q, u, the exact position of the zeros of the function Q, and the 
approximate angles in which the Stokes and anti-Stokes lines emanate from those zeros. 
In the program we use the fact that 



a+dz 



Q{r)dr = -^— {Q{a) + Q'{a)dzf'' - Q{af/' 



(4. 



when the stepsize \dz\ is small enough. Here Q = Q^. The numerical calculation starts at 
the zeros of the function Q. To find the Stokes lines, the program finds the point where 



-e< Re 



2 1 



.3Q'(a) 



{Q{a)+Q\a)dzf/^-Q{af/' 



< e 



and to find the anti-Stokes lines the program finds the points where 



T [2 1 
e < Im < -^^ 



3Q'{a) 



{Q{a) + Q'{a)dzf'-Q{af^ 



< e . 



(4.9) 



(4.10) 



Here e ideally needs to be zero but in the numerical calculations we take it to be a small 
number which generally works well if we take it to be roughly equal to the stepsize \dz\. 
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Figure 5: Stokes (thin) and anti-Stokes (thick) lines produced by numerical calculations 
for pure gravitational axial perturbations, where / = 2, 2M = 1, and q = 0.01. 



As the program calculates the direction of the lines and moves away from the zeros the 
numerical error gets bigger. This is the reason why some times the lines miss each other 
and they move close to each other until the program cannot find any more points as seen 
in Fig. El 

We can use also the numerically generated Stokes and anti-Stokes lines to calculate the 
QNM frequency. First, we use the method of jH] to find the general WKB condition 



g2ir ^ _g2i752 



:i + e 



2J752 ' 



D^hl2 



for the topology shown in Fig. EJ Here 



lab 



h 



Q{r)dr , 



(4.11) 



(4.12) 



and r is given in Eq. ()3.9|) when n = 2. In the Schwarzschild topology, where 712 = 
752 = — TT, Eq. ()4.1H) reduces to the WKB condition we found in Eq. ()3.8p in four 
spacetime dimensions. The phase integrals 712 and 752 can be calculated numerically. For 
the particular case in Fig. Elwe find 712 ~ —3.076 ■ ■ ■ radians and 752 ~ —3.447 ■ ■ ■ radians. 
Therefore we get 

g2ir ^ gS^c^M _ -2.547 ■ • • + il.380 ■ • • . (4.13) 

In other words Re uM ^ ln(2.897 ■ ■ O/Stt. 
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Figure 6: SttMRc cj as a function of — Im uj for pure gravitational axial perturbations, 
where / = 2, 2M = 1, and q = 0.01. 

Similar to Eq. (J4.1H) . we can find a general expression for the WKB condition for the 
topology shown in Fig. El 



.2ire 



l + e''^'^){l + e-''^'^)e 



2i732'\p-2J7i2 „2i732 I 



1 + e 



-2iri 24743 „2i732 



(4.14) 



where Fg and Fj are given in Eqs. (J3.12J) and (J3.13J) respectively when n = 2. To obtain 
the WKB condition (J4.14J) . we have used the fact that 732 = 754 by symmetry. Note that 
the WKB conditions (J4.1H) and (j4.14|) are also valid in higher spacetime dimensions. 
In the limit q/M <C 1, we can write 



Fp ^ -AiTiujM 1 



16M4 



-A-KiuM 



(4.15) 



and 



Ti ^ niuM ( -\- I . 



Therefore, Eq. ()4.14|) can be approximated to be 



^SttujM 



l-(l+e'*^^2)(l + e 



-2i732N„-2J7i2„2i732, 



1 + e 



21743 21732 ^ 



(4.16) 



(4.17) 



We can now numerically calculate the phase integrals 712, 732, and 743 and evaluate the 
real part of oj by using the WKB condition ()4.17|1 . 
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Figure 7: The zeros of the function Q for pure gravitational vector perturbations in 5-d for 
different values of |u;|. 

Fig. iniplots 87rMRe u; as a function of — Im u for pure gravitational axial perturbations, 
where we have taken / = 2, 2M = 1, and q = 0.01. It is interesting to note that at the 
point where the topology changes from R-N to Schwarzschild, the real part of the QNM 
frequency approaches to zero. 



5. Schwarzschild Limit in Higher Spacetime Dimensions 

So far we have concentrated on the metric perturbations in four spacetime dimensions. 
The same type of numerical analysis can readily be done in higher spacetime dimensions 
and one expects similar results. The crucial point is that as the magnitude of the frequency 
u changes, there should be a topology change of the kind which is shown in Figs. El and 
E] This type of topology change happens when two of the zeros of the function Q, which 
are off the real axis initially, coincide on the real axis and then move along the real axis 
staying relatively close to the inner Cauchy horizon. We plot the zeros of the function Q 
for pure gravitational vector perturbations in five spacetime dimensions for different values 
of |u;| in Fig. [7|and give the corresponding results for the real part of the QNM frequencies 
in Fig. ISl respectively. 

We expect this behaviour is generic to all dimensions: in d = n + 2 dimensions there 
are, in the R-N case, two times [n — 1) horizons distributed in pairs along equally spaced 
radial lines in the complex plane. The physical pair is, of course, on the real axis. In the 
Schwarzschild limit Q^ has 2n zeros, while in the R-N limit it has 4n — 2. The difference is 



14 



1.8 
1.6- 

1.4 -" 



ln(3+4cos(27i/5)) 



3 1.2- 
^ 0.81- 



ln(3) 

--■Q ©^ 



0.6- 



0.4- 



0.2 




1- 



10 10 10 10 



9 10 11 12 13 14 

10 10 10 10 10 10 



-Imco 



Figure 8: SvTyuRe a; as a function of — Ini uj for pure gravitational vector perturbations in 
5-ci, where / = 2, 2/i = 1, and 6 = 0.001. 

2{n — 1) zeros, which converge in pairs near the n — 1 inner horizons and then move apart 
at right angles. This behaviour should occur in the range: 



1 < |a;|> 



2„2/(n-l) ^ f ^\ 



4n/(n-l) 



(5.1) 



There will be a critical value of the damping in this range for which the pairs coincide and 
the real frequency will be zero or near zero. In Figs. [7|and|Hl this critical damping happens 
at approximately 1.05 x 10^. 



6. Conclusions 

We have calculated explicitly the QNM frequencies of R-N black holes in the transi- 
tion region of damping that interpolates between the R-N and Schwarzschild limits. Our 
numerical results suggest that the real part of the frequency decreases to zero from its 
ln(3) Schwarzschild value at a critical value of the damping before increasing to the R-N 
value. This behaviour is due to the discontinuous nature of the transition from six evenly 
spaced zeros of the WKB phase to four evenly spaced zeros, with two zeros even closer 
to the origin. In particular the two zeros first meet near the horizon at the critical value 
of the damping and then emerge at right angles. It is at this stage that the Stokes lines 
undergo a discrete topology change and the real part of the frequency goes to zero. We 
have shown that this behaviour is also present in five dimensions, and argued that it is a 
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generic feature of all higher dimensional R-N black holes. 

The physical interpretation of the highly damped QNM frequencies in terms of the 
quantum spectrum of black holes is highly speculative, and seems to be contradicted by a 
myriad of results for multi-horizon, non-asymptotically flat black holes. It should be noted, 
however, that Hod's original conjectures only apply directly to single horizon black holes, 
whose horizons are described by one dimensionful parameter, which could without loss of 
generality be taken as the temperature. The present analysis, in our opinion, seems to 
reverse slightly this negative trend. In effect, we have shown explicitly that for R-N black 
holes with small charge, the infinitely damped QNMs essentially "see" the full structure of 
the solution, including the inner horizon, which affects the value of the frequency in this 
limit. There is nonetheless an intermediate damping for which the QNMs only probe the 
outer horizon, and for this range the "universal" value of ln(3) is reproduced. 

Although our combination of analytic and numerical techniques does not allow us to 
determine whether the real part of the QNM frequency goes exactly to zero, it is tempting 
to speculate that there is a connection between the value of the damping for which this 
happens and the algebraically special frequencies found by Chandrasekhar JH]. In the 
limit of small charge g, one of these values (the one associated with the almost pure 
electromagnetic perturbations) occurs at —i[3l{l + l)M/(4g^)], which is the same order 
of magnitude for which our special frequency occurs. However, our frequency appears 
to be independent of /. Moreover, we have examined in detail only gravitational axial 
perturbations, for which Chandrasekhar's algebraically special frequency gets very small. 
Nonetheless it is important to discover the precise nature of the special frequency that 
occurs in our analysis, and its connection with other algebraically special frequencies. This 
is currently under investigation. 

Hod's physical interpretation of the highly damped QNM frequencies in terms of the 
quantum spectrum of black holes is highly speculative, and seems to be contradicted by a 
myriad of results for multi-horizon, non-asymptotically flat black holes. It should be noted, 
however, that an equally spaced area spectrum for black holes has been derived by several 
authors using a variety of techniques |19j. Moreover, Hod's original conjecture only applies 
directly to single horizon black holes, whose horizons are described by one dimensionful 
parameter, which could without loss of generality be taken as the temperature.^ In any case, 
the present analysis, in our opinion, seems to reverse slightly the current trend, and lends 
support to the possible validity of Hod's conjecture. In effect, we have shown explicitly 
that for R-N black holes with small charge, the infinitely damped QNMs essentially "see" 
the full structure of the solution, including the inner horizon, which affects the value of 
the frequency in this limit. There is nonetheless an intermediate damping for which the 
QNMs only probe the outer horizon, and for this range the "universal" value of ln(3) is 
reproduced. It therefore seems that this issue is worthy of continued investigation. 
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^An intriguing speculation by Makela et. al.[M\ may be relevant in this regard: they propose that for 
multi-horizon black holes it is the total area of all horizons that has an equally spaced spectrum. 
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